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The work is devoted to the investigation of the physical nature 
of H-bonds. The H-bond potential &[j(r,Q) is studied as an ir- 
reducible part of the interaction energy of water molecules. It is 
denned as a difference between the generalized Stillinger-David 
potential and the sum of dispersive and multipole interaction po- 
tentials. The relative contribution of <J?^(r, f2) to the intermolec- 
ular potential does not exceed (10 -r 15)%. 



1. Introduction 

The simplest structure of intermolecular interaction po- 
tential is characteristic of the atomic liquids such as ar- 
gon. Their intermolecular interaction potential <I>(r) is 
the sum of the attractive part $ d i s (f) caused by disper- 
sive forces and $ r0 p(r) that describes the repulsion: 



$(r) = $rep(r) +$di s (V). 



(1) 



In particular, the well-known Lennard- Jones potential 
has such a structure. For molecules having no spherical 
symmetry, the intermolecular potential becomes angular 
dependent [![: 



$(r) ->■ $(r, SI) = $ rcp (r, SI) + $ dis (r, SI) 



(2) 



where SI denotes the set of angles which describe the 
relative orientation of molecules. Such form of the in- 
termolecular potential is characteristic of molecules N2. 
For molecules without center of inversion, it is necessary 
to consider the dipole-dipole interaction and multipole 
interactions of higher order. This circumstance leads to 
the additional term $ M (r,Sl): 



*(r, SI) = $ rcp (r, fi) + $ dis (r, SI) + $ M (r, SI). 



(3) 



The analogous structure of interaction potential is also 
inherent to molecules in water and water-alcohol solu- 
tions if their electronic shells do not overlap. At small 
distances between molecules, the overlapping effect of 
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the electronic shells becomes essential. The correspond- 
ing interaction is usually called the hydrogen bond (H- 
bond). The intermolecular potential is represented in 
the form Q 



$(r,fi) = $ dis (r,tt) + $ H (r,ft), 



(4) 



where $//(r, SI) includes also the repulsion and multipole 
interactions. 

From the qualitative point of view, such a change of 
priorities is not justified, since the analytic continuation 
of multipole contributions to the overlapping region does 
not lead to the effects that can violate the continuity 
condition of the potential. That is why we should rede- 
fine the H-bond potential. In accordance with this, the 
H-bond potential will be defined as 

$(r, SI) = $ rcp (r, SI) + $ dis (r, SI) + $ M (r, SI) + $ ff (r, n), 

(5) 

where $cK S (r, SI) and ^m(t, SI) are the analytic continu- 
ations in the overlapping region of electronic shells. 

One of the most characteristic manifestations of H- 
bonds in liquid water and its vapor is the formation of 
dimers and multimers of higher order. In other words, 
the properties of a dimer give us the direct information 
about H-bonds. This fact points a way how to study 
the intermolecular interaction in water and, rigorously 
speaking, the formation of H-bonds. Let us note the 
main steps of this approach: I) the interaction energy 
of two water molecules is described with the help of the 
most suitable phenomenological model potential that de- 
scribes the ground state of a dimer; 2) the energy ob- 
tained in such a way is compared to that for a water 
dimer calculated with the help of the asymptotic mul- 
tipole expansion; 3) to determine the H-bond potential, 
we construct the difference between the model potential 
and the sum of dispersive and multipole contributions. 
We expect that this difference will be non-zero only in a 
small region near the equilibrium distance between the 
water molecules in a dimer. 
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Fig. 1. Equilibrium configuration of a water dimer 



,.-9-. 




Fig. 2. Definition of angle a 



We suppose that the generalized Stillinger-David po- 
tential (GSD) [3 is the most suitable intermolecular po- 
tential for the description of the interaction between wa- 
ter molecules. This is a soft potential, whose parameters 
can change under the influence of nearby molecules. This 
important circumstance cannot be taken into account 
for almost all phenomenological potentials • In con- 
trast to the original work [9j, the role of the screening 
functions (or overlapping effects) in GSD is taken into 
account more adequately. Moreover, the asymptotic be- 
havior of the original Stillinger-David potential (SD) is 
corrected for large distances between molecules. 

The multipole moments (up to octupole) can be de- 
termined by the quantum-chemistry methods with suf- 
ficient accuracy [10, 11]. Due to this, we are able to 
construct the asymptotic estimation of the interaction 
potential for distances larger than the screening radius. 
At the same time, the analytic continuation of the mul- 
tipole potential to the overlapping region does not lead 
to any serious errors. This allows us to construct the 
above-mentioned difference between the suitable model 
potential and multipole contribution. 

The main goal of this work is the realization of the 
program formulated above for the construction of the 
H-bond potential. 




Fig. 3. Comparative behavior of GSD potential (1), multipole 
potential (2), and SD potential (3) 



2. Ground State of Water Dimer 

In this section, we will represent the results of our study 
of dimers with the help of the GSD potential [3]. In 
order to facilitate the calculations, we will use the hard 
model of water molecules (i.e., the position of hydrogens 
and oxygen remain fixed, as well as their configuration). 
According to [3], such a requirement leads to the error 
of (1.5 -7- 3)%. The ground state of a dimer is identified 
to the minimum of the interparticle potential for two 
molecules presented in Fig. 1. 

The equilibrium distance Tqq between oxygens and 
the angle 9q between the dipole moments of a water 
molecule are determined from the absolute minimum 
of the GSD potential: min $gsd(? : , 9, a = 0). It is 
considered as a function of the dimensionless distance 
f = r oo/ r OH, where toh = 0.97 A is the distance be- 
tween the oxygen and a hydrogen in a water molecule, 
the angle 9 describes some bend of the H-bond, and a is 
related to the rotation around the H-bond (see Fig. 2). 

The comparative behavior of ^gsdj multipole poten- 
tial <1>mp that will be studied below, and $sd from the 
original work [9] is presented in Fig. 3. The equilib- 
rium values of the distance Tqq between the oxygens in 
a dimer, the angle 9q, and the ground-state energy Ed 
(Ed = Ed/kBT mi T m is the melting (crystallization) tem- 
perature for liquid water, fee is the Boltzmann constant) 
for all potentials investigated are presented in Table. In 
addition, the dipole moment Dd of a dimer is included to 
Table as well. It is determined according to the formula 

D d = 2d w cos^(9 + S/2)Ki V2d™ « 2.6D. (6) 

Here, d w denotes the modulus of the water molecule 
dipole moment, S is the equilibrium angle between lines 
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connecting oxygen and hydrogens. The second line in 
Table corresponds to the GSD potential, in which the 
screening length for the function 1 — L[f) (see Q) is 
decreased by A = 0.98 times. The dependence of the 
interaction energy on the rotation angle a around the 
H-bond is especially important (see Fig. 4). 

As is seen, the inequality A$(a) < 1, AQ(a) = 
*^oo Auc^-Hf oo. e a , o) ^ mdicateg the possibility of in _ 

tramolecular rotations of water molecules in a dimer. 
This circumstance should manifest itself in the entropy 
and heat capacity behavior [12]. 

3. Multipole approximation for the interparticle 
potential 

The usage of the multipole interaction potential is jus- 
tified by the following reasons: 1) multipole interaction 
potential of two water molecules is more convenient in 
specific calculations; 2) determination of multipole mo- 
ments within the methods of quantum chemistry is much 
simpler than the construction of approximating func- 
tions; 3) comparison of the attractive part of model po- 
tentials and the electrostatic multipole potential for the 
dimer configuration allows us to control the degree of 
applicability of model potentials. 

The last requirement, as it will be shown below, leads 
to the conclusion that the GSD potential is the most 
suitable for the description of water dimers. 

The multipole potential of the intermolecular interac- 
tion is modeled as 

$mp(0 = $ exp(-fc(r-cr)) +$ M (r), (7) 

where the repulsion part will have the same behavior as 
that for the GSD [3j. Here, tr is the doubled hard-core 
radius of a water molecule, which is identified as the H- 
bond length and $jf(r) is a part of the multipole 
expansion for the interaction energy between two water 
molecules. 



Equilibrium parameters of water dimer 



Dimer parameters 


r oo> 


do, deg 


E d 


d d , D 


GSD 


3.07 


31.1 


-11.36 


2.76 


GSD(A = 0.98) 


3.00 


18.2 


-13.8 


3.02 


<t> M p(5D) 


3.00 


26.8 


-8.73 


3.23 


*mp(5P) 


3.01 


27.9 


-8.19 


3.08 


[13] 


2.98 


60 


-10.69 




[14] 


2.925 


47.5 


-9.11 




[15] 






-10.32 




[16] 


2.976 


57 ± 10 
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Fig. 4. a-dependence of the interaction energy for a dimer at 
r o = 3 A 



We use the multipole expansion $a/(0 up to the 
quadrupole-quadrupole and dipole-octupole terms: 

$Af(r) = + <$>DQ + $QQ + DO + ■■■ ■ (8) 

For the dimer configuration in Fig. 1, the correspond- 
ing terms are 

d 2 

<S>dd = -§— (cos(<5/2 + 6) - 3 cos 8/2 cos 9), 
r oo 

*dq = -^M<£ ) Q$> + &Q$}- 

zr oo 

-i5{4 1} Q$ + 4 a) Q$}). 

$qq = tJ-W^Q^ - 20Q$Q<$ + Q$Q<$), 
4 'oo 

$ do = + 4 2) o&}- 

zr oo 

The analysis of the repulsive part of the multipole po- 
tential shows that the parameters <&o, k, and a take val- 
ues $o = 8.5, k = 5, and a — 2.8. The relative values of 
different contributions to <1>a/ are shown in Fig. 5. Lines 
(3) and (4) testify that the quadrupole-quadrupole and 
dipole-octupole interactions have the same order of mag- 
nitude. 

As is seen from Fig. 3, the multipole part is able to 
correctly describe the interaction energy at distances up 
to 4 A. More exactly, on the interval 4 A < roo < inf, 
the values of multipole potential and GSD coincide with 
high accuracy. At smaller distances between oxygens, 
the overlap of electronic shells begins. The multipole 
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Fig. 5. Role of the terms <$>DDt &DQ, ^QQ, an d &DO m the 
multipole interaction potential (7) of two water molecules in the 
dimer configuration: 1 - &dd, 2- $dq, 3— <3>qq, 4 ~ *BO; 5 — 
■Jrcp, 6 - <J>mp 
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Fig. 6. Partial contributions to the intermolecular potential. 



approach becomes inapplicable. At the same time, the 
phenomenological model potentials are supposed to be 
used also in this region. In particular, the applicability of 
the GSD potential within the overlapping region is jus- 
tified by a suitable selection of the screening functions. 
Unfortunately, the most of other phenomenological po- 
tentials have no necessary compliance. The generalized 
Stillinger-David potential, used in our consideration of 
dimer properties, is quite satisfactory, and it has the 
ability for further modifications. 

4. Hydrogen Bond Potential 

According to the definition of H-bond given in Intro- 
duction, we will consider the difference ^n(r) between 
the generalized Stillinger-David potential <I>gsd and the 
sum of the multipole potential 4>m and the dispersive 
energy H: 

$ ff (r,fi) = $GSD(r,^)-$ r (r,0)-$ dis (r,^)-$ M p(r-,r!). 

(9) 

The behavior of the potentials $gsd(^, ( &dis('", 
and ^Mpf^O) for the configuration of molecules char- 
acteristic of a dimer is presented in Fig. 6. 



0„ 



-2- 




-4- 

Fig. 7. H-bond potential. 

The H-bond potential for the same configuration is 
presented in Fig. 7. 

Hence, &h (r, fl) is a short-range potential that ap- 
pears due to the overlapping of the electronic shells, 
and it has the quantum mechanical nature. It should 
be interpreted as an H-bond potential in water. It 
takes the same order of magnitude as the dispersive 
term and is much smaller than the multipole interaction 
($ dis (l-2)fc B T m , $ M p (7-8)fc B T m , $ H (2~3)k B T m ). 
That is why the contribution of H-bonds to the thermo- 
dynamic potential can be taken into account with the 
help of perturbation theory. This circumstance is con- 
firmed qualitatively by the similarity of thermodynamic 
properties of water and argon on their coexisting curves 

m- 

5. Discussion 

The relatively small depth of the potential well of a 
hydrogen bond leads to the following conclusions: 1) 
contributions of $//(r, fi) to the thermodynamic poten- 
tials and the kinetic coefficients can be calculated with 
the help of perturbation theory; 2) temperature behav- 
ior of the thermodynamic characteristics such as the 
fraction volume or the heat of evaporation is argon- 
like with satisfactory accuracy. The last conclusion is 
confirmed by the results from [18-20]. The H-bond po- 
tential in Fig. 7 corresponds to the equilibrium orien- 
tation of water molecules. There are no restrictions to 
construct $#(r, J7) for all other orientations: all angu- 
lar dependences of $gsd(?", fi) and $mp(»", ^) are well 
known. However, we have to mention that the ther- 
modynamic properties of water are defined by the av- 
eraged intermolecular potential due to the rotation of 
water molecules [20]. It is shown [20] that the potential 
well depth reduces due to the averaging over angular 
variables. This leads to the correction of the argon- like 
dependences of thermodynamic quantities by at most 
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5% [18, 20]. Nevertheless, the contributions of H-bonds 
exhibit themselves in the heat capacity of water (l2| . 
dipole relaxation, and spectral properties [2]. Another 
important circumstance is the necessity to consider the 
influence of the neighbors on the H-bond potential. This 
collective effect will be studied in details in further pub- 
lications. 
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nPHPQUA BOflHEBOrO 3B'5I3Ky Y BOAI 
II.B. MaxAaunyK, M.LI. McuiOMywc, LB. TKmamoK 
P e 3 K) m e 

Po6oTy npiiCBHieHO flOCjiis»ceHHio <pi3HHHo'i npnpo,zni BO^HeBoro 

3B'$I3Ky, HKHH yTBOpiOeTBCH IVlijK MOJieKyjiaMH BOflH. IToTeHniajT 

BOflHeBoro 3B'si3Ky <!> n (r, Q) po3rjiHflaeTbCH hk He3BiflHa HacTima 
eHepri'i B3aeMO,n,ii MOjieKyji bo^h, mo Bii3naHaeTBCH pi3Hnueio Mi}K 
y3arajibHeHHM noTemriajiOM CTijimsscepa Ta /leBi^a i cyMoio no- 
TemriajiiB ^ncnepciiiHoi Ta MyjibTHnojibHo'i B3aeMOflift. IloKa3aHO, 
mp Bi^HOCHa BejinHtma BHecKy <&j](r, fi) He nepeBunrye (10—15)%. 
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